Spin relaxation of electron and hole polarons in ambipolar conjugated polymers.

The charge-transport properties of conjugated polymers have been studied extensively for opto-electronic device applications. Some polymer semiconductors not only support the ambipolar transport of electrons and holes, but do so with comparable carrier mobilities. This opens the possibility of gaining deeper insight into the charge-transport physics of these complex materials via comparison between electron and hole dynamics while keeping other factors, such as polymer microstructure, equal. Here, we use field-induced electron spin resonance spectroscopy to compare the spin relaxation behavior of electron and hole polarons in three ambipolar conjugated polymers. Our experiments show unique relaxation regimes as a function of temperature for electrons and holes, whereby at lower temperatures electrons relax slower than holes, but at higher temperatures, in the so-called spin-shuttling regime, the trend is reversed. On the basis of theoretical simulations, we attribute this to differences in the delocalization of electron and hole wavefunctions and show that spin relaxation in the spin shuttling regimes provides a sensitive probe of the intimate coupling between charge and structural dynamics.

Interlayer Charge Transport in 2D Lead Halide Perovskites from First Principles.

We report on the implementation of a versatile projection-operator diabatization approach to calculate electronic coupling integrals in layered periodic systems. The approach is applied to model charge transport across the saturated organic spacers in two-dimensional (2D) lead halide perovskites. The calculations yield out-of-plane charge transfer rates that decay exponentially with the increasing length of the alkyl chain, range from a few nanoseconds to milliseconds, and are supportive of a hopping mechanism. Most importantly, we show that the charge carriers strongly couple to distortions of the Pb-I framework and that accounting for the associated nonlocal dynamic disorder increases the thermally averaged interlayer rates by a few orders of magnitude compared to the frozen-ion 0 K-optimized structure. Our formalism provides the first comprehensive insight into the role of the organic spacer cation on vertical transport in 2D lead halide perovskites and can be readily extended to functional π-conjugated spacers, where we expect the improved energy alignment with the inorganic layout to speed up the charge transfer between the semiconducting layers.

Transiently delocalized states enhance hole mobility in organic molecular semiconductors.

Evidence shows that charge carriers in organic semiconductors self-localize because of dynamic disorder. Nevertheless, some organic semiconductors feature reduced mobility at increasing temperature, a hallmark for delocalized band transport. Here we present the temperature-dependent mobility in two record-mobility organic semiconductors: dinaphtho[2,3-b:2',3'-f]thieno[3,2-b]-thiophene (DNTT) and its alkylated derivative, C8-DNTT-C8. By combining terahertz photoconductivity measurements with atomistic non-adiabatic molecular dynamics simulations, we show that while both crystals display a power-law decrease of the mobility (µ) with temperature (T) following µ ∝ T -n, the exponent n differs substantially. Modelling reveals that the differences between the two chemically similar semiconductors can be traced to the delocalization of the different states that are thermally accessible by charge carriers, which in turn depends on their specific electronic band structure. The emerging picture is that of holes surfing on a dynamic manifold of vibrationally dressed extended states with a temperature-dependent mobility that provides a sensitive fingerprint for the underlying density of states.

Electronic and Charge Transport Properties in Bridged versus Unbridged Nanohoops: Role of the Nanohoop Size.

In the field of π-conjugated nanohoops, the size of the macrocycle has a strong impact on its structural characteristics, which in turn affect its electronic properties. In this work, we report the first experimental investigations linking the size of a nanohoop to its charge transport properties, a key property in organic electronics. We describe the synthesis and study of the first example of a cyclocarbazole possessing five constituting building units, namely [5]-cyclo-N-butyl-2,7-carbazole, [5]C-Bu-Cbz. By comparison with a shorter analogue, [4]-cyclo-N-butyl-2,7-carbazole, [4]C-Bu-Cbz, we detail the photophysical, electrochemical, morphological and charge transport properties, highlighting the key role played by the hoop size. In particular, we show that the saturated field effect mobility of [5]C-Bu-Cbz is four times higher than that of its smaller analogue [4]C-Bu-Cbz (4.22×10-5 vs 1.04×10-5  cm2 V-1 s-1 ). However, the study of the other organic field-effect transistor characteristics (threshold voltage VTH and subthreshold slope SS) suggest that a small nanohoop is beneficial for good organization of the molecules in thin films, whereas a large one increases the density of structural defects, and hence of traps for the charge carriers. The present findings are of interest for the further development of nanohoops in electronics.

Exciton Dissociation in a Model Organic Interface: Excitonic State-Based Surface Hopping versus Multiconfigurational Time-Dependent Hartree.

Quantum dynamical simulations are essential for a molecular-level understanding of light-induced processes in optoelectronic materials, but they tend to be computationally demanding. We introduce an efficient mixed quantum-classical nonadiabatic molecular dynamics method termed eXcitonic state-based Surface Hopping (X-SH), which propagates the electronic Schrödinger equation in the space of local excitonic and charge-transfer electronic states, coupled to the thermal motion of the nuclear degrees of freedom. The method is applied to exciton decay in a 1D model of a fullerene-oligothiophene junction, and the results are compared to the ones from a fully quantum dynamical treatment at the level of the Multilayer Multiconfigurational Time-Dependent Hartree (ML-MCTDH) approach. Both methods predict that charge-separated states are formed on the 10-100 fs time scale via multiple "hot-exciton dissociation" pathways. The results demonstrate that X-SH is a promising tool advancing the simulation of photoexcited processes from the molecular to the true nanomaterials scale.

Exciton transport in molecular organic semiconductors boosted by transient quantum delocalization.

Designing molecular materials with very large exciton diffusion lengths would remove some of the intrinsic limitations of present-day organic optoelectronic devices. Yet, the nature of excitons in these materials is still not sufficiently well understood. Here we present Frenkel exciton surface hopping, an efficient method to propagate excitons through truly nano-scale materials by solving the time-dependent Schrödinger equation coupled to nuclear motion. We find a clear correlation between diffusion constant and quantum delocalization of the exciton. In materials featuring some of the highest diffusion lengths to date, e.g. the non-fullerene acceptor Y6, the exciton propagates via a transient delocalization mechanism, reminiscent to what was recently proposed for charge transport. Yet, the extent of delocalization is rather modest, even in Y6, and found to be limited by the relatively large exciton reorganization energy. On this basis we chart out a path for rationally improving exciton transport in organic optoelectronic materials.

Charge Transport in Organic Semiconductors: The Perspective from Nonadiabatic Molecular Dynamics.

Organic semiconductors (OSs) are an exciting class of materials that have enabled disruptive technologies in this century including large-area electronics, flexible displays, and inexpensive solar cells. All of these technologies rely on the motion of electrical charges within the material and the diffusivity of these charges critically determines their performance. In this respect, it is remarkable that the nature of the charge transport in these materials has puzzled the community for so many years, even for apparently simple systems such as molecular single crystals: some experiments would better fit an interpretation in terms of a localized particle picture, akin to molecular or biological electron transfer, while others are in better agreement with a wave-like interpretation, more akin to band transport in metals.Exciting recent progress in the theory and simulation of charge carrier transport in OSs has now led to a unified understanding of these disparate findings, and this Account will review one of these tools developed in our laboratory in some detail: direct charge carrier propagation by quantum-classical nonadiabatic molecular dynamics. One finds that even in defect-free crystals the charge carrier can either localize on a single molecule or substantially delocalize over a large number of molecules depending on the relative strength of electronic couplings between the molecules, reorganization, or charge trapping energy of the molecule and thermal fluctuations of electronic couplings and site energies, also known as electron-phonon couplings.Our simulations predict that in molecular OSs exhibiting some of the highest measured charge mobilities to date, the charge carrier forms "flickering" polarons, objects that are delocalized over 10-20 molecules on average and that constantly change their shape and extension under the influence of thermal disorder. The flickering polarons propagate through the OS by short (≈10 fs long) bursts of the wave function that lead to an expansion of the polaron to about twice its size, resulting in spatial displacement, carrier diffusion, charge mobility, and electrical conductivity. Arguably best termed "transient delocalization", this mechanistic scenario is very similar to the one assumed in transient localization theory and supports its assertions. We also review recent applications of our methodology to charge transport in disordered and nanocrystalline samples, which allows us to understand the influence of defects and grain boundaries on the charge propagation.Unfortunately, the energetically favorable packing structures of typical OSs, whether molecular or polymeric, places fundamental constraints on charge mobilities/electronic conductivity compared to inorganic semiconductors, which limits their range of applications. In this Account, we review the design rules that could pave the way for new very high-mobility OS materials and we argue that 2D covalent organic frameworks are one of the most promising candidates to satisfy them.We conclude that our nonadiabatic dynamics method is a powerful approach for predicting charge carrier transport in crystalline and disordered materials. We close with a brief outlook on extensions of the method to exciton transport, dissociation, and recombination. This will bring us a step closer to an understanding of the birth, survival, and annihiliation of charges at interfaces of optoelectronic devices.

Impact of Nanoscale Morphology on Charge Carrier Delocalization and Mobility in an Organic Semiconductor.

A central challenge of organic semiconductor research is to make cheap, disordered materials that exhibit high electrical conductivity. Unfortunately, this endeavor is hampered by the poor fundamental understanding of the relationship between molecular packing structure and charge carrier mobility. Here a novel computational methodology is presented that fills this gap. Using a melt-quench procedure it is shown that amorphous pentacene spontaneously self-assembles to nanocrystalline structures that, at long quench times, form the characteristic herringbone layer of the single crystal. Quantum dynamical simulations of electron hole transport show a clear correlation between the crystallinity of the sample, the quantum delocalization, and the mobility of the charge carrier. Surprisingly, the long-held belief that charge carriers form relatively localized polarons in disordered OS is only valid for fully amorphous structures-for nanocrystalline and crystalline samples, significant charge carrier delocalization over several nanometers occurs that underpins their improved conductivities. The good agreement with experimentally available data makes the presented methodology a robust computational tool for the predictive engineering of disordered organic materials.

Long-Range Interactions Boost Singlet Exciton Diffusion in Nanofibers of π-Extended Polymer Chains.

Raising the distance covered by singlet excitons during their lifetimes to values maximizing light absorption (a few hundred nm) would solve the exciton diffusion bottleneck issue and lift the constraint for fine (∼10 nm) phase segregation in bulk heterojunction organic solar cells. In that context, the recent report of highly ordered conjugated polymer nanofibers featuring singlet exciton diffusion length, LD, in excess of 300 nm is both appealing and intriguing [Jin, X.; et al. Science 2018, 360 (6391), 897-900]. Here, on the basis of nonadiabatic molecular dynamics simulations, we demonstrate that singlet exciton diffusion in poly(3-hexylthiophene) (P3HT) fibers is highly sensitive to the interplay between delocalization along the polymer chains and long-range interactions along the stacks. Remarkably, the diffusion coefficient is predicted to rocket by 3 orders of magnitude when going beyond nearest-neighbor intermolecular interactions in fibers of extended (30-mer) polymer chains and to be resilient to interchain energetic and positional disorders.

Mechanoelectric Response of Single-Crystal Rubrene from Ab Initio Molecular Dynamics.

A robust understanding of the mechanoelectric response of organic semiconductors is crucial for the development of materials for flexible electronics. In particular, the prospect of using external mechanical strain to induce a controlled modulation in the charge mobility of the material is appealing. Here we develop an accurate computational protocol for the prediction of the mechanical strain dependence of charge mobility. Ab initio molecular dynamics simulations with a van der Waals density functional are carried out to quantify the off-diagonal electronic disorder in the system as a function of strain by the explicit calculation of the thermal distributions of electronic coupling matrix elements. The approach is applied to a representative molecular organic semiconductor, single-crystal rubrene. We find that charge mobility along the high-mobility direction a⃗ increases with compressive strain, as one might expect. However, the increase is larger when compressive strain is applied in the perpendicular direction than in the parallel direction with respect to a⃗, in agreement with experimental reports. We show that this seemingly counterintuitive result is a consequence of a significantly greater suppression of electronic coupling fluctuations in the range of 50-150 cm-1, when strain is applied in the perpendicular direction. Thus our study highlights the importance of considering off-diagonal electron-phonon coupling in understanding the mechanoelectric response of organic semiconducting crystals. The computational approach developed here is well suited for the accurate prediction of strain-charge mobility relations and should provide a useful tool for the emerging field of molecular strain engineering.

How to calculate charge mobility in molecular materials from surface hopping non-adiabatic molecular dynamics - beyond the hopping/band paradigm.

Charge transport in high mobility organic semiconductors is in an intermediate regime between small polaron hopping and band transport limits. We have recently shown that surface hopping non-adiabatic molecular dynamics is a powerful method for prediction of charge transport mechanisms in organic materials and for near-quantitative prediction of charge mobilities at room temperature where the effects of nuclear zero-point motion and tunneling are still relatively small [S. Giannini et al., Nat. Commun., 2019, 10, 3843]. Here we assess and critically discuss the extensions to Tully's original method that have led to this success: (i) correction for missing electronic decoherence, (ii) detection of trivial crossings and (iii) removal of decoherence correction-induced spurious charge transfer. If any one of these corrections is not included, the charge mobility diverges with system size, each for different physical reasons. Yet if they are included, convergence with system size, detailed balance and good internal consistency are achieved.

Nonadiabatic dynamics with quantum nuclei: simulating charge transfer with ring polymer surface hopping.

Investigation of many electronic processes in molecules and materials, such as charge and exciton transport, requires a computational framework that incorporates both non-adiabatic electronic effects and nuclear quantum effects, in particular at low temperatures. We have recently developed an efficient semi-empirical fewest switches surface hopping method, denoted fragment orbital-based surface hopping (FOB-SH), that was tailored towards highly efficient simulation of charge transport in molecular materials, yet with nuclei treated classically. In this work, we extend FOB-SH and include nuclear quantum effects by combining it with ring-polymer molecular dynamics (RPMD) in three different flavours: (i) RPSH with bead approximation (RPSH-BA) as suggested in Shushkov et al., J. Chem. Phys., 2012, 137, 22A549, (ii) a modification of (i) denoted RPSH with weighted bead approximation (RPSH-wBA) and (iii) the isomorphic Hamiltonian method of Tao et al., J. Chem. Phys., 2018, 148, 10237 (SH-RP-iso). We present here applications to hole transfer in a molecular dimer model and analyze detailed balance and internal consistency of all three methods and investigate the temperature and driving force dependence of the hole transfer rate. We find that RPSH-BA strongly underestimates and RPSH-wBA overestimates the exact excited state population, while SH-RP-iso gives satisfactory results. We also find that the latter predicts a flattening of the rate vs. driving force dependence in the Marcus inverted regime at low temperature, as often observed experimentally. Overall, our results suggest that FOB-SH combined with SH-RP-iso is a promising method for including zero point motion and tunneling in charge transport simulations in molecular materials and biological systems.

Quantum localization and delocalization of charge carriers in organic semiconducting crystals.

Charge carrier transport in organic semiconductors is at the heart of many revolutionary technologies ranging from organic transistors, light-emitting diodes, flexible displays and photovoltaic cells. Yet, the nature of charge carriers and their transport mechanism in these materials is still unclear. Here we show that by solving the time-dependent electronic Schrödinger equation coupled to nuclear motion for eight organic molecular crystals, the excess charge carrier forms a polaron delocalized over up to 10-20 molecules in the most conductive crystals. The polaron propagates through the crystal by diffusive jumps over several lattice spacings at a time during which it expands more than twice its size. Computed values for polaron size and charge mobility are in excellent agreement with experimental estimates and correlate very well with the recently proposed transient localization theory.

Crossover from Hopping to Band-Like Charge Transport in an Organic Semiconductor Model: Atomistic Nonadiabatic Molecular Dynamics Simulation.

The mechanism of charge transport (CT) in a 1D atomistic model of an organic semiconductor is investigated using surface hopping nonadiabatic molecular dynamics. The simulations benefit from a newly implemented state tracking algorithm that accounts for trivial surface crossings and from a projection algorithm that removes decoherence correction-induced artificial long-range charge transfer. The CT mechanism changes from slow hopping of a fully localized charge to fast diffusion of a polaron delocalized over several molecules as electronic coupling between the molecules exceeds the critical threshold V ≥ λ/2 (λ is the reorganization energy). With increasing temperature, the polaron becomes more localized and the mobility exhibits a "band-like" power law decay due to increased site energy and electronic coupling fluctuations (local and nonlocal electron-phonon coupling). Thus, reducing both types of electron-phonon coupling while retaining high mean electronic couplings should be part of the strategy toward discovery of new organics with high room-temperature mobilities.

Detailed balance, internal consistency, and energy conservation in fragment orbital-based surface hopping.

We have recently introduced an efficient semi-empirical non-adiabatic molecular dynamics method for the simulation of charge transfer/transport in molecules and molecular materials, denoted fragment orbital-based surface hopping (FOB-SH) [J. Spencer et al., J. Chem. Phys. 145, 064102 (2016)]. In this method, the charge carrier wavefunction is expanded in a set of charge localized, diabatic electronic states and propagated in the time-dependent potential due to classical nuclear motion. Here we derive and implement an exact expression for the non-adiabatic coupling vectors between the adiabatic electronic states in terms of nuclear gradients of the diabatic electronic states. With the non-adiabatic coupling vectors (NACVs) available, we investigate how different flavours of fewest switches surface hopping affect detailed balance, internal consistency, and total energy conservation for electron hole transfer in a molecular dimer with two electronic states. We find that FOB-SH satisfies detailed balance across a wide range of diabatic electronic coupling strengths provided that the velocities are adjusted along the direction of the NACV to satisfy total energy conservation upon a surface hop. This criterion produces the right fraction of energy-forbidden (frustrated) hops, which is essential for correct population of excited states, especially when diabatic couplings are on the order of the thermal energy or larger, as in organic semiconductors and DNA. Furthermore, we find that FOB-SH is internally consistent, that is, the electronic surface population matches the average quantum amplitudes, but only in the limit of small diabatic couplings. For large diabatic couplings, inconsistencies are observed as the decrease in excited state population due to frustrated hops is not matched by a corresponding decrease in quantum amplitudes. The derivation provided here for the NACV should be generally applicable to any electronic structure approach where the electronic Hamiltonian is constructed in a diabatic electronic state basis.

Electron and excitation energy transfers in covalently linked donor-acceptor dyads: mechanisms and dynamics revealed using quantum chemistry.

Photoinduced electron transfer (ET), hole transfer (HT), charge recombination (CR) and energy transfer (EET) are fundamental mechanisms, which occur in both natural and artificial light harvesting systems. Here, we present a computational strategy which determines ET, HT, CR and EET rates in a consistent way and merges them in a kinetic model to reproduce the net excited state dynamics. The effects of the solvent are included in all steps of the calculations making the present strategy a useful tool for a rational design of charge and energy transfer processes in complex systems. An application to covalently linked zinc and free-base porphyrin-naphthalenediimide dyads is presented. For each of the two systems, ultrafast optical spectroscopy experiments have shown a specific photophysics with different processes taking place simultaneously. The model reveals that such a diversity is mainly due to the different relative stability of the charge-separated state, while the electronic couplings for charge and energy transfer processes are quite similar in the two dyads.